Differential Photoelectron Holography: 
A New Approach for Three-Dimensional Atomic Imaging 
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We propose differential holography as a method to overcome the long-standing forward-scattering 
problem in photoelectron holography and related techniques for the three-dimensional imaging of 
atoms. Atomic images reconstructed from experimental and theoretical Cu 3p holograms from 
Cu(OOl) demonstrate that this method suppresses strong forward-scattering effects so as to yield 
more accurate three-dimensional images of side- and back-scattering atoms. 

PACS numbers; 61.14.-x, 42.40.-i 



Holography P] is a method of recording both the am- 
plitudes and phases of waves scattered by an object il- 
luminated with coherent radiation, and using this infor- 
mation to directly construct a three-dimensional image 
of the object. Szoke first suggested that coherent out- 
going waves from atomically-localized sources of photo- 
electrons, fluorescent x-rays, and 7-rays could be used 
to achieve atomic-scale holography. This idea was ini- 
tially demonstrated theoretically for the case of photo- 
electrons by Barton 3] , and then extended into a multi- 
energy format by Barton and Terminello and by Tong 
and co-workers '4|. By now several experimental ap- 
proaches to such atomic- resolution holo^a^hy have been 
demonstrated, including photoelectrons 5, 6, 7, 8], Auger 
electrons j^J , Kikuchi electrons ^3 j diffuse-scattered low- 
energy electrons 0| , fluorescent x-rays in either a direct 
mode |l2j or a multi-energy inverse mode 7-rays 
[l^. and bremsstrahlung x-rays 15]. 

Among these methods, photoelectron holography (PH) 
has the advantages of being capable of studying the lo- 
cal atomic structure around each type of emitter without 
requiring long-range order and of distinguishing emitters 
through core-level binding-energy shifts 0. Photoelec- 
tron holograms also show strong modulations of up to 
±50%, so such effects are easily measurable. However, 
PH can suffer from serious image aberrations due to the 
strength of electron scattering. The atomic scattering 
factor / is a highly anisotropic function of scattering an- 
gle, and can depend strongly on electron kinetic energy 
Ek- In particular, as Ek increases above a few hundred 
eV, / becomes more and more significant in the forward 
direction, resulting in a strong forward-scattering (FS) 
peak [1^ that can induce image aberrations. Beyond 
this, PH also can suffer from multiple-scattering (MS) 
effects due to the scattering strength. 

Various reconstruction algorithms and measurement 



methods 0, IH S ^| have been proposed to correct for 
the anisotropic / and MS effects, some of which can be 
summarized via 



U{v) = 



j Wx (k) cxp [-ikr + ik ■ r] d^k 



(1) 



where U is the image intensity at position r, x is the 
normalized 3D hologram, and the function or operator 
W permits describing the difference between algorithms, 
with W=l in the original multi-energy formulations 0|. 
One alternative algorithm sets W = {k,9^) so 
as to divide out the anisotropic /, where 6'Jf is the an- 
gle between r and k. In another algorithm 7] based on 
the more ideal electron back scattering (BS), a window 
function for W that limits the integral in Eq. to be 
in a small cone of k around — r is chosen to emphasize 
the imaging of BS atoms. Although successful in several 
applications 0, difficult to apply this small- 

cone method to many systems where the imaging of FS 
or even side-scattering (SS) atoms is important, such as 
epitaxial films and buried interfaces. In fact, imaging of 
"bulk" atoms surrounded by FS and BS atoms via PH 
has proven to be especially difficult [cf. Figs. 7-9 in ref. 
0], with most successful applications being to emitters 
in the first few layers near a surface. 

To overcome such FS effects, we propose in this Let- 
ter "differential holography" . By simply replacing x in 
Eq. Q by its /c-derivative (i.e. W — d/dk) or more 
conveniently by a numerical difference between two x's 
at different energies {Sx = + Sk) — ef- 
fects can be greatly suppressed. We have applied this 
method to multi-energy holograms for Cu 3p emission 
from Cu(OOl), and show that this provides images that 
are improved over prior work in several respects. 

To avoid confusion with other methods in PH, we also 
note that "derivative" PH has been proposed and used 
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successfully by Chiang and co-workers However, the 
purpose here is to eliminate uncertainties in I due to the 
variation of experimental conditions by first taking loga- 
rithmic derivatives [dl /dk]/I that are then reintegrated 
into "self-normalized" intensities; thus, it is still finally x 
that is used in Eq.l^l. 

The principle of differential photoelectron holography 
(DPH) is as follows. We consider the single-scattering 
expression of x for an emitter-scatterer pair spaced by a 
vector r 

X(k)= ^ 

« MM^cos [kr (1 - cose^) + if (fc, 9^)] , (2) 

where /q is the intensity that would be observed without 
atomic scattering, and ip is the scattering phase. If 5k 
is sufficiently small so that 5 | / 1 / 1 / 1 ^ 1 , where 5\f\ 
is the change in |/|, the difference of two holograms at 
k± — fc ± (5fc/2 can be written in a similar form to Eq. 
as: 

« -ilMl sin [kr (1 - cos6i^) + if {k, 6^)] , (3) 

where direction k is defined by angles 6 and 0, the 
"effective" scattering amplitude is defined as |/gff| = 
2 I/I sin [Skr (l — cos0jf)/2 -I- Sip/2\ , and (p is the average 
of if's at k±. In the PS region where 6^ — > 0, |/gff| is 
thus very small, approaching zero in the limit of Sip 0. 
If Sk is also small, l/gffj is proportional to r; thus, DPH 
not only suppresses the FS effects, but also enhances the 
imaging of distant atoms. In Pig. ^ |/| and |/eff| are 
plotted as a function of 9^ for Cu-Cu nearest neighbors 
(r=2.56 A). For k=A.6 A^^ and A-\ j/g^j is 

significant only in the region of 9^ >~ 90°. Therefore, 
the imaging of SS and BS atoms is expected, while it 
will be difhcult for this case to image FS atoms. On the 
other hand, for k=8.8 A~^ and a larger fractional ^A;=1.0 
A~^, |/gff| is significant not only in the BS region but 
also in the range of ^30°-90°. Since near-neighbor 
PS diffraction fringes extend out beyond 30° 0, 0| , we 
might expect the latter choice to also permit imaging FS 
atoms. In this way, the relative sensitivity of DPH to SS 
and FS atoms can be "tuned" by selecting the range and 
step width of k scans. Finally, we note that the suppres- 
sion of MS effects by means of a transform over a volume 
in k space is well known in normal multi-energy PH |J| 
and this suppression will be equally present in DPH. If 
anything, the inherent elimination of strong FS effects in 
DPH should lead to even better MS suppression. 

To demonstrate DPH experimentally, photoelectron 
holograms from Cu(OOl) were measured at beamline 7.0 
of the Advanced Light Source at the Lawrence Berke- 
ley National Laboratory. Photoelectron spectra for Cu 



3p emission were collected at 25 energies over fc=4.5- 
9.3 A-i {Ek=77-330 eV) with a constant step of 5k=0.2 
A^^ (SEk =7-14 eV), along 65 different directions over a 
symmetry-reduced 1/8 of the total solid angle above the 
specimen, and with a polar angle range from ^=0° (sur- 
face normal) to 70°. The photoelectron intensity I{k,9,(f)) 
was fitted by low-order polynomials to obtain the smooth 
background intensity [a |2i| : 

/o (fc, 9) = (ao + aifc 02^^) (60 + hi cos 9 + b2 cos 36*) . 

(4) 

Three kinds of x were obtained from this fitting: xa 
by fitting the second factor of Eg. (4) to a scanned-angle 
pattern Ik [9, (f>) at each fixed fc |a|j Xs by fitting the first 
factor to a scanned-energy curve I^ (k) at each fixed di- 
rection k and xc by fitting both factors to the full 
data set of I (fc, 9, cj)) at one time, with the last expected 
to be the most accurate from an a priori point of view 
[SI]. The fc-differences from xc were also used for DPH 
in what we will term Method D (i.e., xd — ^xc)- Since 
low-frequency fringes due to FS events in I^ (fc) are au- 
tomatically removed in Method B |0, the resulting Iq 
inherently deviates from the true Iq defined as the in- 
tensity without scattering, especially in the FS direction. 
In addition, since Ik {9, <p) and I^ (fc) are independently 
normalized without considering the continuity of x iii 
the whole sampled k space in Methods A and B, they 
could degrade holographic fringes in /j^ {k) and Ik {9, (/<), 
respectively. By contrast. Method C takes into account 
the continuity of x over the whole data set, but the FS 
peaks remain in xcJ however, they should be eliminated 
in Xd ■ The original transform of Eq. ^ was used for all 
four data sets; but to avoid the abrupt truncation of the 
integral in Eq. (0), W was taken to be the product of 
a Gaussian function of k and a Hanning function cos^ 9, 
with an additional multiplication by r to make atoms at 
larger distances more visible. 

Figure El shows atomic images reconstructed from XA- 
Xd in the vertical (100) plane of Cu(OOl). In Methods A 
and C, only elongated features related to FS effects from 
atoms of types 6 and 7 are observed above Zc = —0.5 
A (the arbitrary location of a change in image multipli- 
cation). This is consistent with a previous PH study of 
W(llO) 0], in which Method C was used. By contrast, it 
has been reported P that FS atoms of type 5 have been 
imaged via Method A from Cu 3p holograms for Cu(OOl) 
obtained at 9 energies. Even if possible differences in the 
two sets of experimental data are taken into account, it 
is difficult to conclude from our results that the images 
of these FS atoms can be resolved from strong artifacts 
via Method A. Below Zc, several peaks near the BS po- 
sitions 1-3 are observable for A and C among various 
strong artifacts, but only with the help of higher image 
amplifications of 46 and 29, respectively. 

In Methods B and D, image intensities are stronger 
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FIG. 1: Comparison of the usual scattering amplitude |/| and 
the effective scattering amplitude of differential holography 
|/eff|' calculated for Cu-Cu nearest neighbors (r=2.56 A) as 
a function of scattering angle Sj^ for two different sets of k 
and &k in taking the differential of x- i.^) ^ ~ 4-6 (81 
eV), &k = 0.2 A"^ (7 eV) and (b) k = 8.8 A'^ (295 eV), l>k 
= 1.0 A~^ (67 eV). The final strong forward-scattering data 
points of I/I at the right of panel (b) are truncated. 



in the BS region, with the relative image amplification 
factors being reversed in sense and smaller at x5 com- 
pared to A and C. In Method D, a strong, somewhat 
elongated peak is observed at the FS position 6, with 
weaker features that appear to be associated with atoms 
7 also present in the corners of the image. In both B 
and D, two strong peaks are observed at the SS positions 
4 above Zc and five peaks are observed at the BS posi- 
tions 1-3 below Zc- However, the most intense features 
in Method B are the artifacts between the two nearest 
BS atoms of type 3. In Method D, by contrast, the five 
strongest peaks below Zc are of roughly equal intensity 
and correspond reasonably well to the near-neighbor BS 
atoms. Therefore, we find Method D to be the most ro- 
bust for imaging both SS and BS atoms (as well as to 
some degree also FS atoms 6), even if there are shifts in 
position of approximately 0.1 A for type 1, 0.6 A for 2, 
and 0.3 A for 3. Such peak shifts relative to the true 
atomic positions, as observed in all methods, can be at- 
tributed to the present neglect of corrections for both the 
scattering phase and the inner potential. 



(a) Notiiial-Jot'f?, <.:i)-Ex. (b) Nomial-J"of1fcLJ-Ex. 
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FIG. 2: Atomic images in the vertical (100) plane of Cu (001) 
reconstructed from Cu 3p holograms obtained by Methods 
A-D, as described in the text. The emitter and scatterer 
positions are indicated by squares and circles, respectively, 
and various near-neighbor atoms are numbered. Image in- 
tensities above or below Zc — tO.5 A have been rescaled by 
the factor shown in each panel, with this factor being de- 
termined so as to make the maximum intensities above and 
below Zc equal. Experimental images (a) Image from Method 
A: normal holography with Iq determined by fitting the an- 
gular variation of Eq. (4) at each k value, (b) Image from 
Method B: normal holography with Jo determined by fitting 
the k variation along each direction, (c) Image from Method 
C: normal holography with Iq determined by fitting both the 
angular and k variation, (d) Image from Method D: differen- 
tial holography, with Jo as in (c). Theoretical images: (e) As 
(c) but theoretical. 

For comparison with experiment, we have also per- 
formed MS simulations of /(k), using a cluster method 
fully described elsewhere |2l|. The theoretical Jq was 
obtained simply as the square of the zeroth-order wave 
function without scattering. Images reconstructed from 
the theoretical x and Sx via Methods C and D are shown 
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in Figs, life) and (f) and can be compared with Figs. [Sfc) 
and (d), respectively. The main features in Figs. HJc) 
and (d) are well reproduced by our simulations, although 
the artifacts between the atoms 3 are much stronger in 
experiment for C, and the relative intensity in the re- 
gion of FS atom 6 is stronger in experiment for D. Even 
though the ideal x was used for image reconstruction, 
no atomically-resolved SS or FS peaks are observable in 
Fig. I3c). Therefore, the corresponding artifacts in Fig. 
EIc) are not purely due to the uncertainties in the exper- 
imental data and any errors in the /q subtraction, but 
must have their origin in the MS effects and basic imag- 
ing algorithm. On the other hand. Fig. |2If) exhibits 
well-resolved peaks at the SS and BS positions. Since 
there are no artifacts below Zc in Fig. Eff), the artifacts 
in Fig. Eld) are by contrast considered to be purely due 
to the experimental noise and other non-idealities in the 
data analysis. 

We have also generated full three-dimensional atomic 
images from the experimental data via X-D; although 
length limitations prevent showing these here. In these 
images, we find in addition to the atoms 1-4 and 6 in Fig. 
13 two other types of near-neighbor BS and SS atoms lo- 
cated in the vertical (110) plane (denoted types 2' and 
4' and situated in the same horizontal layers as 2 and 
4, respectively). All of these atoms are reasonably well 
reconstructed, with only a few, such as 2, being signifi- 
cantly shifted in position, but most within a few tenths 
of an A of the correct positions in all directions. The 
overall positional errors for all of the atoms compared to 
the known Cu lattice can be summarized as (radial lo- 
cation shift in xy)/(vertical location shift in z), and are: 
0.0 A/0.1 A for atoms 1, 0.6 A/0.1 A for 2, 0.3 A/0.1 A 
for 2', 0.2 A/0.1 A for 3, 0.1 A/0.0 A for 4, 0.3 A/0.0 A 
for 4', and 0.0 A/0.4 A for 6. As a further indication of 
the overall image quality obtained by DPH, the reader is 
referred to an animated comparison of 3D images for the 
four approaches of Figs. 2(a)-(d), in which DPH is alone 
in imaging approximately 15 near-neighbor atoms [2^ . 

Finally, we compare DPH with a very recently intro- 
duced approach for PH termed near-node holography 
[2^ , in which FS effects are suppressed by using a special 
experimental geometry with electron exit nearly perpen- 
dicular to light polarization. Although this technique is 
promising, DPH has the advantages that it does not re- 
quire a special experimental geometry or s-subshell-like 
form for the photoelectric cross section, that it seems to 
yield images of as good or better quality [i^, 113 and 
that it can be used in other types of holography in which 
polarization cannot be varied. 

In summary, we have demonstrated differential pho- 
toelectron holography (DPH) as a powerful method for 
overcoming the FS problem in PH and enhancing image 
quality for any kind of system in which FS can arise, as 
for example, bulk emission and buried interfaces. This 
method should also be helpful in other types of electron 



holography in which energy can be stepped in a controlled 
way (e.g. Kikuchi jlOj] or LEED 'll'l holography). The re- 
constructed images for Cu 3p/Cu(001) demonstrate that 
DPH is successful in suppressing the FS effects so as to 
image SS, BS, and to some degree also FS, atoms with 
accuracies of 0.1-0.6 A. 
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